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ABSTRACT 

It is shown that the metal abundance in the solar interior is constrained from the 
current solar neutrino experiment, even if one allows neutrino conversion in the sun due 
to neutrino oscillation. The result shows that the metal abundance in the interior should 
be within the range 0.4-1.4 times that in the surface, supporting the idea that the sun 
formed by a contraction of a gas cloud with an almost homogeneous composition. 
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1. Introduction 

Whether the metal abundance measured in the 
stellar surface represents correctly the value in the 
stellar core is a nontrivial question. In computations 
of the evolution of stars it is generally assumed that 
the initial metal abundance of the core is identical to 
the present day abundance of the surface, for which 
spectroscopic information is available for heavy ele¬ 
ments. This agrees with the simple picture that stars 
formed by a contraction of the molecular cloud of the 
uniform composition. 

In principle, however, the metallicities of the core 
and the surface can be different, unless the whole stars 
are fully convective like in low mass (M < OAMq) 
main sequence stars (e.g., Iben 1967). An extreme 
example is Jupiter, where the core consists mostly of 
silicates, while the surface is dominated by hydrogen 
and helium with small admixture of heavier elements. 
It may well be that the origin of Jupiter like stars is 
completely different from that of ordinary stars (e.g., 
Podolak, Hubbard & Pollack 1993): nevertheless it 
would not be absurd to imagine that rocks are taken 
into the core at the time of gas contraction, and the 
core has metallicity higher than in the surface since 
the birth of the star, for which the radiative transport 
dominates. There is also an effect that heavy metals 
tend to sink by diffusion towards the core, while the 
envelope becomes hydrogen rich. A calculation shows 
that such an effect do exist in the sun as evidenced 
from helioseismology, although the amount is small 
(Bahcall et al. 1997). 

The main reason that one assumes the homoge¬ 
neous composition throughout the star is basically 
due to lack of information of the metal abundance 
in the stellar interiors. In this article we consider 
the problem that the sun might offer a possibility to 
study the metal abundance in the deep interior or in 
the core through the solar neutrino observations that 
probe the core region of the sun. This problem, how¬ 
ever, is not very obvious due to the famous solar neu¬ 
trino problem that the observed neutrino flux is by 
a factor 2-4 times smaller than is predicted with the 
standard model of the sun (Bahcall 1989). The most 
elegant and widely accepted solution to this problem 
is that electron neutrinos emitted in the nuclear re¬ 
actions of the sun are converted by neutrino oscilla¬ 
tion into muon or tau neutrinos that are sterile in the 
nuclear detector or have a much smaller cross section 
for scattering off electrons, as advocated by Mikheyev 


and Smirnov (1986a; refereed to as the MSW effect). 

This solution at first glance appears to make the 
things so flexible that almost any amount of the neu¬ 
trino flux before the oscillation effect is experimen¬ 
tally allowed, if one tunes the neutrino mass and mix¬ 
ing parameters in some appropriate way. Indeed, for¬ 
getting about all knowledge on the nuclear reaction 
cross sections, even the case that almost 100 % of 
the solar energy is generated by the CNO cycle is 
not excluded solely by the solar neutrino experiment 
(Bahcall, Fukugita & Krastev 1996). 

This, however, is too extreme. If we adopt the 
knowledge of nuclear reactions within the range al¬ 
lowed by the current experiment, the freedom is not 
that large. For instance, the amount of ®B must be in 
the range between 1 /3-2 times the value the standard 
solar model predicts, and the core temperature can be 
determined to within 5% of the standard solar model 
value (Hata & Langacker 1997; hereafter HL97) in or¬ 
der to satisfy the current solar neutrino experiments, 
whatsoever the neutrino parameters one takes. In 
this paper we study to what extent the current solar 
neutrino experiments constrain the metal abundance 
in the sun allowing for MSW neutrino oscillation (or 
neutrino conversion). 

The increase of the metal abundance obviously pro¬ 
motes the CNO cycle. It also increases opacity that 
modifies the core temperature. Therefore, it increases 
the relative importance of the neutrino flux from the 
CNO cycle significantly. This increment of the CNO 
neutrino flux must be cancelled by increasing impor¬ 
tance of the suppression factor coming from the neu¬ 
trino oscillation to keep the consistency with the solar 
neutrino experiments. In so far as this works, a larger 
metal abundance in the solar interior is allowed. If the 
metal abundance is increased more, however, we no 
longer have solutions that satisfy all the solar neutrino 
experiments in a consistent way. This is the logic that 
we are going to explore in this study. 

We remark that metal abundance may not be so 
tightly constrained from helioseismology information. 
Helioseismology is sensitive to the sound velocity, 
the change of which reads approximately Acs/cs ~ 
1/2(AT/T — A/r/^) with T and /x temperature and 
mean molecular weight, = 2A-|-0.75y-|-0.58Z(^^C)-|- 
.... A cancellation takes place between temperature 
and mean molecular weight. Although the present 
best solar model is known to give Cg as accurate as 
0.2%, and the 50% change in Z might be detected 
if it changes with satisfying AZ = —AX, the pres- 
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ence of V complicates the situation. Hence, it is not 
obvious how strong constraint can be derived on the 
metallicity independently of X and Y. 

2 . Calculation 

We take the standard solar model of Bahcall and 
collaborators as our fiducial (Bahcall & Ulrich 1988; 
Bahcall & Pinsonneault 1992; 1995, hereafter BP95), 
and consider a small departure from their best model. 
The energy of the sun is generated from the pp chain 
(98%) and the CNO cycle (2%). Neutrinos are pro¬ 
duced in pp dev, pep —> dv, e+'^Be^'^Li+v and 
®B^®Be-|-e-|-i^ from the pp chain, and in beta decay of 
13 n^ isq and i7p jn tPe CNO cycle. The experimen¬ 
tal information comes from the long-running Homes- 
take experiment with ^^Cl, which is sensitive to both 
high energy neutrinos (®B neutrinos) and intermedi¬ 
ate energy neutrinos (^Be, pep and CNO neutrinos) 
(Cleveland et al. 1997), a water Cerenkov experiment 
at Kamiokande and Superkamiokande measuring only 
for high energy neutrinos (Fukuda et al. 1996; Tot- 
suka et al. 1997), and gallium experiment, Gallex 
(Hampel et al. 1996) and Sage (Abdurashitov et al. 
1996) that are very sensitive to low energy neutrinos 
{pp neutrinos). The problem is that the detection 
rate is smaller than predicted, by factors, 3.7 ± 0.6, 
2.6±0.5 and 2.0±0.2, respectively. Furthermore, this 
specific energy dependent suppression pattern makes 
the explanation of the problem by modifying the solar 
model highly unlikely, leaving the neutrino oscillation 
explanation as the most attractive possibility (Bahcall 
& Bethe 1990; Fukugita & Yanagida 1991). 

Indeed, this energy dependent suppression is very 
naturally understood by considering the conversion of 
electron neutrinos into other types of neutrinos inside 
the sun by neutrino oscillation (Mikheyev & Smirnov 
1986a). For the neutrino flux given by the standard 
solar model, the current solar neutrino experiments 
allow two small parameter regions that are located in 
the two parameter space, neutrino mass-square dif¬ 
ference Ato^ = {i = p or r) and intrinsic 

mixing angle {0) between the two relevant neutrinos: 
one is called the small angle solution, in which in¬ 
termediate energy neutrinos are suppressed, and the 
other the large angle solution, for which the suppres¬ 
sion of neutrino fluxes is almost energy independent. 
The most up to date calculations are found in Bahcall 
& Krastev (1996) and in HL97. 

Additional information can be obtained from an 


upper limit on the possible flux variation between the 
day and night time (day-night effect) (Mikheyev & 
Smirnov 1986b; see Bahcall & Krastev 1997 and refer¬ 
ences therein). For some specific neutrino parameter 
range the converted muon neutrinos are changed back 
to electron neutrinos during the propagation through 
the Earth, causing an increase of the neutrino capture 
rate in night in the Kamiokande and Superkamioka¬ 
nde detectors. The absence of this effect down to 2% 
level (Fukuda et al. 1996) excludes a substantial size 
of parameter regions of our interest. 

We repeat the neutrino propagation calculation al¬ 
lowing for a variation in the metal abundance of the 
sun. We use the scaling law of Bahcall & Ulrich 
(1988), which has given the explicit metallicity de¬ 
pendence for each component of the neutrino flux. 
The range of the model explicitly studied covers about 
±50% around the value of metallicity determined for 
the solar surface. The dependence outside this range 
is a simple extrapolation with power law. Although 
this calculation is rather old, we expect the gross 
metallicity dependence does not differ from what one 
could obtain from the more modern BP95 calcula¬ 
tion. We impose a luminosity constraint so that the 
luminosity that would change due to a change in the 
metal abundance is renormalized to the today’s lu¬ 
minosity of the sun. Namely, we study the model at 
fixed luminosity. 

The table of Bahcall & Ulrich shows that the most 
sensitive to metal abundance are indeed the CNO 
neutrino fluxes: the power 7 of the flux (j) oc {Zj 
is 1.86 for the ^^N neutrino, 2.03 for the neutrino 
and 2.09 for the ^^F neutrino. This high power is 
caused by a multiplicative effect due to the increase 
of the abundance of catalysing ^^C and the increase 
of opacity that makes the core temperature higher. In 
spite of its sharp temperature dependence, the effect 
on the ®B neutrino flux is smaller (7 = 1.27) than for 
the CNO neutrinos. 

As for the Hducial flux, we use the value of the 
BP95 calculation with metal diffusion effect taken 
into account: (('(®B) = ( 6.6 ± 1 . 1 ) x 10 ® cm“^s“^, 
9.3± 1.3 SNU for captures with ^”^01 and 137±8 SNU 
for captures with ^^Ga. We take 2.55±0.14 ± 0.14 
SNU for the Homestake experiment (Cleveland et 
al. 1997), (2.80 ± 0.19 ± 0.33) x 10®cm-2s-i for 
Kamiokande (Fukuda et al. 1996) and (2.51^® :[3 ± 
0.18) X 10®cm“^s“^ for Superkamiokande [combined 
(2.586 ± 0.195) X 10®cm-2s-i] (Totsuka 1997), and 
69.5 ± 6.7 SNU for combined Sage and Gallex exper- 


3 



iments. The absence of the day-night effect is also 
imposed on our data analysis. In our actual calcu¬ 
lation we use the night-time flux data divided into 
5 bins according to the cosine of angle from the sun 
(Fukuda et al. 1996). 

The data are then fitted with the three free param¬ 
eters, Am^, sin^ 20 and ZjX, and calculate a likeli¬ 
hood function taking account of both experimental er¬ 
rors of neutrino reaction rates and those arising from 
uncertainties of solar models as given by BP95, in the 
same way as done in HL97. 

The resulting probability distribution is displayed 
in Fig. 1 taking Z as a parameter. The range allowed 
at 95% confidence level (CL) is 

0.4 < Z/Zsurface < 1.4 (1) 

where ^surface = 0.0175 (Grevesse & Noels 1993) us¬ 
ing the standard solar model value X = 0.71. The 
region outside this range is excluded even if we as¬ 
sume the flux reduction due to neutrino oscillation. 
Namely, the metal abundance in the solar interior 
cannot be much different from that in the surface. 
The allowed range is, of course, much larger than the 
change of the metal abundance induced by the dif¬ 
fusion effect, which is about 15% in Z (BP95). We 
remark that the range given in eq. (1) is the range 
where explicit solar model studies are made by Bah- 
call and Ulrich (1988) and the behaviour regarding 
the variation of Z is well studied; so hindsight we need 
not to use power law of the Z dependence out to the 
range where its behaviour is not well established. 

We also show a figure (Fig. 2) of allowed region for 
the two other parameters, Am^ and sin^ 20, overlaid 
on the corresponding figure with the standard case 
(i.e., Z = ^surface allowing for errors for the value at 
the surface 6.1%). The contour is a parameter range 
corresponding to 95% CL. Compared with the con¬ 
tour of the standard model case, the allowed region is 
elongated horizontally for the small angle solution, 
or vertically for the large angle solution. Most of 
the elongated parts (a part in the left hand side of 
the small angle solution, and that in the upper part 
of the large angle solution) correspond to the case 
with metallicity lower than the default value. For 
larger metallicity, the change required for the neutrino 
parameters is quite small to accommodate increased 
neutrino fluxes. Nevertheless, there is a sharp cut off 
against the increase of metallicity, beyond which ap¬ 
propriate neutrino parameters do not exist to make 
the flux consistent with the three solar neutrino ex- 
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Fig. 1.— Probability distribution as a function of 
metal abundance of the interior of the sun. Likelihood 
is normalized to the maximum height being unity (so 
that 95% CL is at about 0.18). 

periments. In our calculation the constraint from the 
absence of the day night effect serves to squeeze the 
high metallicity end in the large angle solution, but 
it plays little role for the small angle solution. 

3. Discussion 

We have shown that the solar neutrino experiments 
put a strong constraint on the metal abundance in 
the interior of the sun, even if we allow the neutrino 
oscillation due to the matter effect. Namely, the pos¬ 
sibility of neutrino oscillation does not lend us much 
freedom to increase metallicity inside the sun. The 
metallicity in the interior of the sun should not be 
larger than the surface value by more than 40%, or 
0.15 dex in [Fe/H] for the initial value. This is a good 
news to the people working for stellar evolution cal¬ 
culations, since we expect that the sun is not a special 
case, but the same probably applies to more general 
cases, justifying the standard assumption that stars 
formed by a contraction of a homogeneous gas sphere 
(Hayashi 1966). 

On the other hand, the errors in the current neu¬ 
trino experiments and uncertainties in nuclear reac¬ 
tion cross sections still allow the possibility that the 
metal abundance in the solar interior is slightly larger 
(smaller) than in the surface. Accepting this uncer¬ 
tainty, energy generation from the CNO cycle may 
be as uncertain as 0.4% to 4% of total energy gen- 
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Fig. 2.— Allowed neutrino parameters (95 % 
CL) when metallicity is treated as a free parame¬ 
ter (shaded regions). The thick solid contours are 
the region allowed for the standard solar model {Z = 
^surface) at 95 % CL. The contours indicate the sup¬ 
pression factor for the neutrino flux in steps of 
0.1 unit, as discussed in the discussion section below. 


eration, which is compared to 1.8% for the standard 
value (see Fig. 3). The or neutrino fluxes 
may vary from 1x10® to 12xl0®cm“^s“^ within this 
range. 

One might think that the or neutrino flux 
(the end point energies 1.199 or 1.732 MeV, respec¬ 
tively) itself can give us a useful indicator for the 
metallicity in the solar interior due to its sensitivity 
to the carbon abundance. One may prepare a de¬ 
tector with a detection threshold set just above the 
energy of ^Be neutrinos (0.862 MeV). Once a new de¬ 
tector, which can measure ^Be neutrinos, e.g., with 
a liquid scintillator measuring for ve ve, is con¬ 
structed [e.g., Borexino (Arpesella et al. 1992)], this 
is not difficult, since the CNO neutrino flux is 100 
times higher than ®B neutrino flux in the standard 
solar model. The oscillation effect, however, makes 
the situation somewhat subtle. For the allowed re¬ 
gions of neutrino parameters, we expect the MSW 
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Fig. 3.— Fraction of the solar energy borne by the 
CNO chain as a function of the metal abundance in 
the solar interior. 

suppression factor for the ^®N neutrino flux to be 0.1- 
0.5 (0.4-0.7) for the small (large) angle solution for a 
given metallicity. Unless this large uncertainty aris¬ 
ing from neutrino oscillation is reduced, the flux 
does not give us useful information on the metallicity. 

We find that the most effective way to reduce this 
uncertainty is to measure gallium capture rate as pre¬ 
cise as possible. If the error attached to gallium cap¬ 
ture rate is reduced, the region shrinks in the ver¬ 
tical (horizontal) direction for the small (large) an¬ 
gle solution, i.e., the error in the suppression factor 
is reduced. A ®B neutrino flux measurement via the 
neutral current interaction planned at Sudbury, when 
combined with a better gallium experiment, serves 
to further reduce the error in the suppression fac¬ 
tor for a given metallicity (although this experiment 
alone is not very effective for our purpose). For in¬ 
stance, if one achieves an error of the gallium ex¬ 
periment as small as ±2 SNU around the current 
best value and measure the ®B neutrino neutral cur¬ 
rent reaction to a 10% accuracy, we would obtain 
0.85 < Z/^surface < 1-2 at 95% CL, if the small angle 
solution is right, or 0.9 < ^/^surface < 1-1 if the large 
angle solution is correct (see Fig. 4). 

We have argued that solar neutrino experiments 
have given unique information for the metals in the 
solar interior, which is not accessible by other means. 
Assuming that the sun is not special, this removes 
our worry that metal abundance in the stellar inte- 
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Fig. 4.— and ®B neutrino fluxes as a function of 
the metal abundance in the solar interior. Solid and 
dashed curves indicate those fluxes without modifica¬ 
tions due to neutrino conversion (the flux of ®B neu¬ 
trinos is multiplied by a factor of 100). The shaded 
regions are those expected for flux after neutrino 
conversion for the small angle (SA) and large angle 
(LA) solution. We assume that the error of gallium 
capture rate was reduced to ±2 SNU and the ®B neu¬ 
trino neutral current reaction was measured with an 
error of 10 %. 

rior might be different from that in the surface. For 
example, this can be an issue in deriving metallic- 
ity dependence of the Cepheid period-luminosity re¬ 
lation, for which the metal effect on the colour, which 
is affected by surface metal abundance, cancels to a 
large degree against that on the luminosity, which is 
affected by metal abundance in the core of the sun 
(Stothers 1988). If metallicity would be different in 
these two places, the large cancellation no longer takes 
place, resulting in a much larger metallicity depen¬ 
dence in the Cepheid period-luminosity relation. The 
present analysis, if our result applies to star forma¬ 
tion from molecular clouds in general, implies that 
the metallicity difference, at least, is not a likely pos¬ 
sibility to account for a large metallicity dependence 
of the Cepheid period luminosity relation suggested 
recently (Gould 1994; Sasselov et al. 1997; Sekiguchi 
& Fukugita 1997). We have also discussed that one 
could reduce the error in Z/Zgurface to a level of 10- 


20 % with further improvement in some specific solar 
neutrino experiments. At this level, the effect of metal 
diffusion could be seen with the solar neutrino experi¬ 
ments, opening a possibility to carry out a cross check 
with the result from helioseismology. 
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